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Abstract. We investigate the possibility to acquire information on the generalized parton distribution E 
and, through a model for E, also on the it-quark total angular momentum J u by studying deeply virtual 
Compton scattering and hard exclusive p° electroproduction on a transversely polarized hydrogen target 
at Hermes. It is found that a change in J u from zero to 0.4 corresponds to a 4a (2a) difference in the 
calculated transverse target-spin asymmetry in deeply virtual Compton scattering (p° electroproduction), 
where a is the total experimental uncertainty. 
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1 Introduction 

Over more than two decades, inclusive and semi-inclusive 
charged lepton scattering has been used as a powerful tool 
to successfully study the longitudinal momentum struc- 
ture of the nucleon, which was parameterized in terms of 
parton distribution functions (PDFs). Hard exclusive re- 
actions can be described in the theoretical framework of 
generalized parton distributions (GPDs) y2J|2EK|- Their 
application became apparent after it had been shown [S] 
that measurements of the second moment of the sum of 
the 'unpolarized' GPDs H and E open, for the first time, 
access to the total angular momentum of partons in the 
nucleon: 




Pi 



P2 



Fig. 1. In the parton picture, GPDs describe correlations 
between two partons with different longitudinal momenta at 
given Q 2 and i, where t — (pi — P2) 2 also contains transverse 
degrees of freedom. 



J a (Q 2 ) = -Ihn 



H a {x,tt,Q 2 ) + E a (x^,t,Q 2 ) 



dx. 



(1) 

In this relation H a (x,£ l ,t,Q 2 ) and E a (x,£ i ,t,Q 2 ) denote 
parton spin non-flip and spin flip GPDs (a = u,d,s), re- 
spectively 1 . GPDs depend on the fractions x and £ of lon- 
gitudinal momentum of the proton carried by the parton 
and 011 f = (pi — P2) 2 , the square of the 4-momentum 
transfer between initial and final protons (see. Fig.Q. As 
ordinary PDFs, also GPDs are subject to QCD evolution. 
Their Q 2 dependence has been perturbatively calculated 
up to next-to-leading order in a s [B] and is omitted in the 
notations throughout the paper. 



1 Throughout this paper the GPD definitions of a recent re- 
view are used. 



Recently, a simultaneous description of the transverse 
spatial and the longitudinal momentum structure of the 
nucleon was shown to be an appealing interpretation of 
GPDs pilOlllllliijj . The concept of GPDs covers several 
types of processes, ranging from inclusive deeply inelas- 
tic lepton scattering to hard exclusive Compton scatter- 
ing and meson production. Measurements of GPDs are 
expected to shed light especially on the hitherto theoreti- 
cally uncharted territory of long-range ('soft') phenomena 
where parton-parton correlations are known to play an 
important role. 

First steps towards the extraction of the GPD H have 
already been performed by scattering leptons off unpo- 
larized protons through measurements of either cross sec- 
tions |13U14| , or cross section asymmetries with respect to 
beam charge or beam spin | 16U17| . Future measure- 
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ments of the transverse target-spin asymmetry (TTSA) in 
hard exclusive electroproduction of a real photon (deeply 
virtual Compton scattering, DVCS) or a vector meson of- 
fer the possibility to acquire information on the spin-flip 
GPD E. The most promising experiments to access it are 
those running at intermediate energy, where the spin-flip 
amplitude is expected to be sizable, while at higher en- 
ergies it is suppressed due to s-channel helicity conserva- 
tion. Thus at present a realistic program may be envis- 
aged for Hermes, Clas and possibly Compass. In this 
paper the prospects are discussed for Hermes measure- 
ments of TTSAs in DVCS and p° electroproduction, and 
in particular their sensitivity to the u-quark total angular 
momentum. 



2 Modeling Generalized Parton Distributions 

GPDs are most commonly parameterized using an ansatz 
based on double distributions [T^lFTT?j complemented with 
the D-term [201 ■ Factorizing out the i-dependence, the 
non-forward GPDs can be related to the ordinary PDFs 
and the proton elastic form factors. In this framework \'21\ . 
the spin non-flip GPD H is given by 



l-{l + Kp )t/4m 2 H q . g (x,Q 
l-t/4m 2 (l-t/0.71) 2 



where k p =1.793 is the proton anomalous magnetic mo- 
ment and m is the proton mass. The neutron Dirac form 
factor is neglected compared to the one of the proton. 

For quarks, the i-independent part of the GPDs H q is 
written as 

H q (x,t) = H? D (x,t) + 9(£-\x\)D q (£j , (3) 

where D q {^j is the D-term, and H q D is the part of the 
GPD that is obtained from the double distribution (DD) 
F q - 

1 1-1/31 

H^ D {x,0 = J dp J da 6(x-p-at) F q {p,a) . (4) 

-1 -l + l/3| 

For the double distributions the suggestion of Ref. JS| is 
used, 

F q (f3 ) a) = h((3,a)q(f3), (5) 
where the profile function is given by |19j : 



HP, a) 



r{2b + 2) [(1- 



2 2b + 1 r 2 (b+ 1) (1-1/31)26+1 



(6) 



For p > 0, q(f3) — q va i{P) + q(P) is the ordinary quark 
density for the flavor q. The negative (3 range corresponds 
to the antiquark density: q(—/3) — —q(J3). The parameter 
b characterizes to what extent the GPD depends on the 
skewness £. In the limit b — > oo the GPD is independent 



on £, i.e., H(x,£,) = q(x). Note that & is a free parameter 
for valence quarks (b va i) or sea quarks (b sea ) and thus can 
be used as a fit parameter in the extraction of GPDs from 
hard electroproduction data \22\ . 

For gluons, the t- independent part of the GPD H g is 
directly given by the double distribution, 

H g {x,i) = H g DD {x,i) = 

1 1-1/31 

dp J daS(x- P- a£) pF g (P, a) (7) 

-1 -1 + 1/51 

with the same form of the profile function in the double 
distribution 



F g (p,a)=HP,a)g(p). 



(8) 



The t-dependence for gluons is taken to be the same as 
that for quarks. 

The factorized ansatz (j2J is the simplest way of model- 
ing GPDs. However, experimental studies of elastic diffrac- 
tive processes indicate that the ^-dependence of the cross 
section is entangled with its dependence on the photon- 
nucleon invariant mass |2.'ij . Recent evidence comes from 
lattice QCD calculations P^IT25] and phenomenological 
considerations [2323 The non-factorized ansatz can be 
based on soft Regge-type parameterizations. In this case, 
the ^-dependence is not factorized out and not controlled 
by a form factor as in Eq. (J2J). Instead, it is kept in Eqs. Q), 
(@J and |0. The ^-dependence of double distributions is 
then modeled as |21|: 



F q , g {P,a,t) 



1 



(9) 



which is referred to as Regge ansatz in the following. Here 
a 1 is the slope of the Regge trajectory, a' q = 0.8 GeV -2 
for quarks and a' g = 0.25 GeV -2 for gluons. 

In the factorized ansatz the spin-flip quark GPDs E q 
are given by [2*T] : 



E q (x,t,t) = 



E q (x,Q 
(1 - i/0-71) 2 



(10) 



In the Regge ansatz the t-dependence is modeled in anal- 
ogy to Eq. Q. 

The t-independent part is parameterized using the dou- 
ble distribution ansatz: 



E q {x,t)=E q JU {x,t)-6(fi-\x\)D q 



(11) 



Note that the D-term has the same size, but the opposite 
sign in Eqs. Ijll|l and ©• Therefore, it drops out when 
calculating J q according to Eq. (JIJ . 

The double distribution has a form analogous to the 
spin-nonflip case: 



1 1-1/31 

E q D (x,0=Jdp J da5(x-p-a£) K q (P,a) (12) 
-i -i+l/3| 
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with: 



K q (fi,a) = h(fi,a)e q (p). 



(13) 



The spin- flip parton densities e q (x) can not be ex- 
tracted from deep-inelastic scattering (DIS) data, unlike 
the case of spin non-flip ones. Based on the chiral quark 
soliton model [2] , the spin- flip density is taken as a sum of 
valence and sea quarks contributions. Since in this model 
the sea part was found to be very narrowly peaked around 
x = 0, the whole density is written as: 



z q {x) = A q q va i(x) + B q S(x) 



(14) 



In this expression, the shape of the valence quark part is 
given by that of the spin non-flip density. The coefficients 
A q and B q are constrained by the total angular momen- 
tum sum rule (Tj^) and the normalization condition 



/ 



dx e q (x) 



(15) 



where K q is the anomalous magnetic moment of quarks of 
flavor q (k u = 2n p + n n = 1.67, Kd = k p + 2n n = —2.03). 
The constraints yield: 



A„ 



B,, 



2 Jo 



M, 



(2) 



M, 



(2) 
9»«l 



1 _ 2J U - M ( u ] 
2 Ku ^(2) 



( 2) 

Here M y q > and M y q J al 

are the parton momentum contribu- 
tions to the proton momentum: 



B d = Hd - 



r( 2 ) 



2 J d - M 



(2) 



M. 



(2) 



(16) 
(17) 
(18) 



M< qli = I xq va i(x)dx, 



M< 2 > = 

q 



x [q va i(x) + 2q(x)] dx. 

(19) 

In the given scenario the total angular momenta carried 
by u- and ci-quarks, J u and Jd, enter directly as free pa- 
rameters in the parameterization of the spin-flip GPD 
E q (x,^,t). Hence the parameterization (|l"4"}l can be used 
to investigate the sensitivity of hard electroproduction ob- 
servables to variations in J u and Jd- 

As to the gluons, there exists no hint how the spin-flip 
GPD E g could be described. There is an expectation that 
Eg is not large compared to E u and Ed US]. Hence, for 
simplicity throughout the present study E g is neglected 
("passive" gluons, i.e. E g — 0). 

As an example, Fig. [5] shows the i-independent part 
of various GPDs at f = 0.1, based on the MRST98 gS| 
parameterization of PDFs at Q 2 =4 GeV 2 . Using instead 
CTEQ6L PDFs [301 as input, the results for u(d) quark 
GPDs are changed by less than 3%(10%); the GPD H g 
is up to 40% larger at x = 0. Because of u-quark domi- 
nance in electroproduction, uncertainties originating from 




Fig. 2. t- independent part of quark and gluon GPDs at 
Q 2 =4 GeV 2 , £=0.1 (MRST98 PDFs are used). 



d-quark PDFs can be safely neglected. Since gluons are ab- 
sent in leading-order DVCS, uncertainties resulting from 
gluon PDFs are of little influence for DVCS asymmetries 
and have been found to lead to a fractional change of up to 
15% for the p° asymmetries. For the following calculations 
the MRST98 PDF set is taken. 



3 Sensitivity of DVCS to the w-quark Total 
Angular Momentum 

3.1 Cross Section and Asymmetries 

The 5-fold cross section for the process e(fc)+p(pi)— ►e(fc')- 
p(j>2) + 7(92) is given by: 



da 



a em X BV 



dx B dQ 2 dtd<t>d<t> s 1&t: 2 Q 2 ^/T 



Ax 2 B m 2 /Q 2 



where Q 2 = — q\ is the negative squared 4- momentum of 
the virtual photon, xb = Q 2 / (2pi • qi) is the Bjorken vari- 
able, t — (pi — P2) 2 , y — {pi ■ <Zi)/(pi ' k), T denotes the 
photon production amplitude and e is the electron charge. 
Since the DVCS and Bethe-Heitler (BH) processes have an 
identical final state, in which the photon is radiated either 
from a parton or from a lepton, respectively, T is given 
by the coherent sum of the BH amplitude Tbh and the 
DVCS amplitude T DVC s- 



\T\ 



IT, 



BH 



Tdvcs\ 2 = \T BH \ 2 + \T DV cs\ 2 +X, (21) 



in which 



I = TX M T, 



BH 2 DVCS 



XbhTdvgs 



(22) 



describes the interference between both processes. 

The coordinate system is defined in the target rest 
frame, as explained in Fig. The theoretical formulae 
used below refer to the target being transversely polar- 
ized w.r.t. the virtual photon direction, while in the ex- 
periment the target polarization is transverse w.r.t. the 
incident lepton direction. At Hermes kinematics, these 
two directions are approximately parallel and the small 
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Fig. 3. Kinematics and azimuthal angles of photon electro- 
production in the target rest frame. The z-direction is chosen 
along the three-momentum of the virtual photon qi. The lep- 
ton three-momenta k and k' form the lepton scattering plane, 
while the three-momenta of virtual and real photons q\ and q2 
define the production plane. The azimuthal angle of the pro- 
duction plane with respect to the scattering plane, around the 
virtual photon direction, is denoted as <f). Correspondingly, rf>s 
denotes the azimuthal angle of the target polarization vector 
with respect to the lepton scattering plane. In this frame the 
target polarization vector is given as S± = (cos 4>s, sin (f>s, 0). 
The definitions conform with the Trento conventions |31|. 



longitudinal component (< 10%) of the target polariza- 
tion along the virtual photon direction can be neglected. 
Thus the reasonable approximation 



da = da unp + daTP 



(23) 



is used, where da unp (daTp) denotes the cross section for 
the unpolarized (transversely polarized) component. 

Since in the kinematic region of the Hermes exper- 
iment the DVCS cross section is typically much smaller 
than the BH cross section [32], the contribution of the 
DVCS term to the total cross section is neglected in the 
following. The contributions of the BH term for an unpo- 
larized beam are: 



\T BH I 

| unp | 



\n-BH\ 
\ I TP 



x 2 B y 2 (l + Ax 2 B m 2 IQ 2 ) 2 tPx{c\>)P 2 ^) 

Kunp + C?,unp COS <\> + c|f np COS 2<j>] , (24) 

0. 



The full expressions for the BH propagators P\ (</)), P2{4>) 
and for the Fourier coefficients cf^ np can be found in 
Ref. E2] 2 . 

The leading twist and leading order a s contributions of 
the DVCS-BH interference term to the total cross section 
can be written as: 



_ ±e b 



T 



±e b 



TP 



xlyHP x {ct>)P 2 {ct>) 



(r 1 -4- r 1 

\ 0,tmp * °1 ,imp 

/ (x Bl y,Q 2 ) x 



cos< 



(25) 



IiuMn sin (0 — <j>s) cos <p + ImM$ cos {(p — <f>s) sin <f> 



The azimuthal angles defined in this work are different from 



those used in Ref. 



(f> = 7T 



'33 



and 



Here + (— ) stands for a negatively (positively) charged 
lepton beam and f \x B ,y, Q 2 ) is a kinematic pre-factor 
independent of azimuthal angles. The full expressions for 
c\ unp can be found in Eqs. (53-56) of Ref. [231 - Mn and 
M$ are certain linear combinations of the Compton form 
factors H, £, H and £ , which are convolutions of the re- 
spective twist-2 GPDs H, E, H and E with the hard- 
scattering kernels as defined in Eq. (9) of Ref. 

The full expressions for Mm and Ms can be found 
in Eq. (71) in Ref. [S| or in Eq. (60) in Ref. E2|. Since 
£ ~ xb/{2 — xb) is small in a wide range of experimentally 
relevant kinematics, terms with pre-factor ^ or x B can be 
neglected, except for the GPD E because the pion pole 
contribution to E scales like £ , so that Mm and Ms can 
be approximated as: 



M N 



M, 



t 



AM 2 
t 

'AM 2 



F 2 H -FiS], 
F 2 H - F^£ 



(26) 



Here F\ and F 2 are the Dirac and Pauli form factors of 
the proton, respectively. 

In order to constrain the GPDs involved in Eq. I|26l) . 
the transverse polarization component of the interference 
term, Ttp, has to be singled out. This can be accomplished 
by forming the transverse (T) target-spin asymmetry with 
unpolarized (U) beam: 



A 



r tf, - (j) S ) = da ^ ~ ^ ~ da ^ ~ ^-s + 7r ) (27) 
da{<j) — 4>s) + da((f> — <j>s + n) 

~ Ap T ^ s) cos • sin(0 - S ) cos 



+ ^cos(0-0 s ) S in^ cos( 



4>s)sm<t>. (28) 



As | Tj^p | and 1 unp are independent on — 0s, they do 
not appear in the numerators of Eq. I|27|l . Since their dom- 
inant contribution to the denominator in Eq. I|27|) is given 

by cgg np , the two amplitudes of the TTSA, A s ^ 

■ >sm( ^, can be approximated as: 



(<t>-4>s) 



and A 



UT 



A 



UT 



±f (xb,v,Q 2 



ImM 



N 



.cos (<t>—<t>s) sin <j> 
l UT 



±f(x B ,y,Q 2 )- 1 -^. (29) 



"0,unp 



Note that the approximations used in this section are 
for illustrative purposes only and are not used in the nu- 
merical calculations described below. 



3.2 Expected Value of TTSA and Projected Statistical 
Uncertainty 

Since in the DVCS process the gluons enter only in NLO 
in as, their contributions to cross section and TTSA are 
neglected. For the quarks, it can be seen from Eq. Q26JI 
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that, besides the GPDs H and E which have been dis- 
cussed in Sect. there are two other GPDs, H and E, 
involved in the TTSA for DVCS. Since they are not the 
main interest of this paper, in the calculations below they 
are always included and kept unchanged. In their model 
description, the forward limit of the GPD H is fixed by 
the quark helicity distributions Aq(x,fi 2 ), while the GPD 
E is evaluated from the pion pole which only provides a 
real part to M s in Eq. 

At present, there exists a code [35] designed to cal- 
culate observables in the exclusive reaction ep — > ejr/. It 
has been used (see App. EJ) to evaluate the TTSA arising 
from the DVCS-BH interference. The TTSA is calculated 
at the average kinematic values per bin in ig, Q 2 and t 
taken from a measurement of the beam-spin asymmetry 
in DVCS at Hermes (see Tab-CJ). 

The statistical error of an asymmetry is independent 
on its size if the asymmetry itself is small. For a single 
beam (target) spin asymmetry it is obtained as: 



1 

N 



1 



p2 

beam{tar get) 



(30) 



where N is the total number of events that is proportional 



to the integrated luminosity, and P bt 



^(target' 



is the beam 



(target) polarization. The following projection is based on 
a future Hermes data set of 8 million DIS events to be 
taken with an unpolarized positron beam and a trans- 
versely polarized hydrogen target. Using the known sta- 
tistical errors of the beam-spin asymmetry measurement 
at Hermes on an unpolarized hydrogen target (7 million 
DIS events, Pbeam — 50%) [32], the projected statistical 
error for the TTSA is obtained. 



Table 1. Average kinematic values for Q 2 , xb, —t bins and 
statistical errors, taken from a measurement of the beam-spin 
asymmetry at Hermes |22| . 



Q 1 bin (GeV-') 


1.00-1.50 


1.50-2.30 


2.30-3.50 


3.50-6.00 


6.00-10.0 


(Q') (GeV^) 


1.2 


1.8 


2.8 


4.4 


7.1 


{xb) 


0.06 


0.08 


0.10 


0.15 


0.24 


(-*) (GeV^) 


0.07 


0.09 


0.12 


0.17 


0.24 


stat. &A a J ffl > 


0.053 


0.050 


0.061 


0.070 


0.163 


xq bin 


0.03-0.07 


0.07-0.10 


0.10-0.15 


0.15-0.20 


0.20-0.35 




1.4 


2.2 


3.1 


4.5 


6.1 


{xb) 


0.05 


0.08 


0.12 


0.17 


0.24 


(-t) (GeV 2 ) 


0.08 


0.10 


0.12 


0.17 


0.22 


stat. SA^ll^ 


0.048 


0.053 


0.060 


0.099 


0.145 


-t bin (GeV^) 


0.00-0.06 


0.06-0.14 


0.14-0.30 


0.30-0.50 


0.50-0.70 




2.0 


2.5 


3.0 


3.6 


3.9 


{x B ) 


0.08 


0.10 


0.11 


0.12 


0.12 


(-t) (GeV^) 


0.03 


0.09 


0.20 


0.37 


0.57 


stat. SA^l* 


0.041 


0.052 


0.066 


0.126 


0.263 



-0s) COS0 



and A 



cos {<fi—(fis) sin tf> 



The projections for <mu i f ; 

are calculated for different values of the total angular 
momentum J u . Since the contributions of u-quark and 
d-quark are proportional to the corresponding squared 
charge, the d-quark contribution is suppressed and hence 
in the calculations a fixed value is used for Jd- The latter 



was chosen to be Jd = 0, inspired by the results of recent 
lattice calculations (see e.g. Ref. [HE]). Using both Regge 
and factorized ansatze, the asymmetries are calculated for 
the four possible cases setting the profile parameters b va i 
and b sea to either one or infinity. Comparing all sets of 
projections to each other, the amplitudes of the TTSA 
appear to be sensitive only to the change in b sea from one 
to infinity. The resulting differences are small and can be 
seen by comparing Figs. 0] and |5j where the amplitudes 
are shown in dependence on Q 2 , xb and —t together with 
the projected statistical errors. In order to study the con- 
tributions of the GPDs H, H and E alone, calculations 
are done for E = as well. 




-t(GeV') 



-t(GeV') 



Fig. 4. Expected DVCS TTSA amplitudes A% (4 " 4 ' s) cos and 
^cos(</> ^> s )sm<j> j n Regge ansatz for b va i = 1, b sea = oo, 
J u = 0.4 (0.2, 0.0), J d = 0.0. E = denotes zero effective con- 
tribution from the quark GPDs E q . The calculations are done 
at the average kinematic values as listed in Tab. Projected 
statistical errors are shown. 



As expected from Eqs. (|26|l and l|29(l . variations in 
the parameter settings for the GPD E become manifest 



A 



sin {(p—^>s) cos (j> 
UT 



while A 



cos {<j)—(ps) sin <f> 
UT 



shows only minor 



modifications. The latter are apparent only in the kine- 
matic regime of large xb or correspondingly large Q 2 since 
the contribution of the GPDs E q to Ms is suppressed by 
xb and thus has been neglected in Eq. I|26(l . Within these 



model calculations A 



sin {4>—<t>s) cos 4> 
UT 



turns out to be sizable 



even when the calculation is done for E q — 0. Thus a solid 
knowledge about the GPD H u is needed in order to con- 
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0.2 



i 



-0.2 



4 6 8 

Q 2 (GeV 2 ) 



-0.4 
-0.6 



: , 4 : 






: J u =o.4 

J>0.2 ... 


J =0 



4 6 8 

Q 2 (GeV 2 ) 




0.4 0.6 

-t(GeV 2 ) 



t(GeV') 



Fig. 5. Expected DVCS TTSA amplitudes A^ T 

.cos (<j>-4> s ) sin <j> . 



sin (<^) — (fig) cos 



and 



A UT in the Regge ansatz for b v 



al — 1, bsea — 1; 

J u = 0.4 (0.2, 0.0), J d = 0.0. E = denotes zero effective con- 
tribution from the quark GPDs E q . The calculations are done 
at the average kinematic values as listed in Tab. Projected 
statistical errors are shown. 



strain J u . It has been shown \AT\ that the model parame- 
ters for the GPD H u , in particular the size of the profile 
parameters b va i and b sea , can be well constrained by the 
envisaged Hermes DVCS measurements until 2007, using 
an unpolarized hydrogen target. Since in addition the pro- 
file parameters are assumed to be the same for the GPD 
E u , the only remaining free parameter is J u . Hence the 
projected measurement of jf^.'t>-'t> s ) cos 't> j ias a c ' e ar po- 
tential to constrain J u , as can be seen from the left panels 
of Figs. H and 

The discriminative power of the envisaged TTSA mea- 
surement can be enhanced by combining all data into one 
point, since all considered models show the same kine- 
matic dependences. The corresponding statistical power 
of a Hermes data set based on 8 million DIS events is 
shown in Fig. |SJ for b sea equal to one or infinity, and for 
three different values of the total u-quark angular momen- 
tum J u plus the special case E q = 0. 

It appears that for b sea — 1 (b sea — oo) the ampli- 
tude ranges between values of -0.17 and -0.27 (-0.19 and 
-0.29) when J u ranges between 0.4 and 0. The projected 
statistical error for these integrated TTSA amplitudes is 
0.017. Extrapolating the knowledge on the systematic un- 
certainty from the analysis of 2000 Hermes data [22], its 
size can be expected to not exceed the statistical error, 



IV -0.2 
< 



-0.3 





J u =°-4 - 








J u =0.2 j 




j„=o i 
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Fig. 6. Expected DVCS TTSA amplitudes ^° w_< * s) cos 
with b va i — 1 and b sea = oo (left panel) or b sea = I (right 
panel), J u — 0.4 (0.2, 0.0), Jd = 0.0 in the Regge ansatz at the 
average kinematics of the full measurement. E = denotes 
zero effective contribution from the GPDs E q . The projected 
statistical error for 8 million DIS events is shown. The system- 
atic error is expected to not exceed the statistical one. 



such that a total experimental uncertainty below 0.025 
appears as a realistic estimate. Altogether, the difference 
in the size of the TTSA due to a change of J u between zero 
and 0.4 corresponds to a 4<r effect, where a denotes the 
total experimental uncertainty. Thus, based on the GPD 
model used it can be expected that the upcoming DVCS 
results from Hermes 3 will provide a constraint on the size 
of J u . 



4 Sensitivity of Elastic p° Electroproduction 
to the u-quark Total Angular Momentum 

Also a measurement of the TTSA in elastic vector meson 
electroproduction can be a source of information about 
the spin-flip generalized parton distribution E. An esti- 
mate for the asymmetry was obtained in Ref. |21| . using 
the factorized model of GPDs described in Sect. 2 with- 
out inclusion of gluons. The scope of this section is to also 
include the Regge ansatz, to check the assumption that 
the gluon contribution to the p° electroproduction cross 
section is small, and to eventually calculate the size of the 
TTSA at Hermes kinematics. The issue is raised since 
in contrast to DVCS, in vector meson elastic electropro- 
duction gluons enter at the same order of a s as quarks, 
namely at order a s to the power one. Hence this channel 
appears as one of the rare cases where gluon GPDs may 
be accessed through Hermes data. 



4.1 Cross Section and Gluonic Contribution 

It was shown j88| that the leading twist contribution to ex- 
clusive electroproduction of vector mesons requires both 
the virtual photon and the vector meson to be longitu- 
dinal, i.e. transversely polarized. Therefore the present 



3 The recent switch of the accelerator (HERA) to an electron 
beam will require to also perform the above calculations for 
the negative beam charge. However, the sensitivity to J u of 
the combined electron and positron measurements is expected 
to be similar to the one calculated here for a positron beam 
only. 
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calculations cover only the longitudinal part of the cross 
section. 

The cross section of the reaction ^* L p — > p°p' is given 

by EH 



1 



d(JL _ 

dt ~ 8rmr(W 2 - m 2 )\ qi \ 



(\Ta\ 2 + \T b 



(31) 



where is the momentum of the virtual photon in the 
center of mass system of this photon and the initial proton, 
while W is their invariant mass. The spin-flip amplitude 
reads 

l 

_ 2\/27ra s . , . ,. . x f , $(z) ,„„, 

9Q J z 

o 

~ -iAirea s ^-^- [ dz^—^- 7 
9Q J z 
o 



and the spin- flip one is 4 : 
\[2"K(X 



Mz) 



T R = e^-^Bu(p 2 )a^n^AMPi) / (33) 
9<ym J z 

o 

to 9 C 



Here n = (1, 0, 0, — l)/(V2(pi + 2?2) + ) is a light-like vector 
along the z-axis, Z\ = j»2 — Pi is the 4-momentum trans- 
fer (A 2 = t). The modulus of its transverse component 
is given by \At\ = yj—t(l — £ 2 ) — 4£ 2 m 2 . The p°-meson 
wave function is taken in the form 



= 6z(l - z)f p 



(34) 



with / p =0.216 GeV and z being the meson longitudinal 
momentum fraction carried by a parton. The complex fac- 
tors A and B are given by: 



A = 



1 

J= J (e u H u (x,£,t) - e d H d (x,£,t) 



o(e u - e d ) 



{ i . + ; . } 



(35) 



B = J (e u E u (x,£,t) - e d E d (x,£,t) 



;(e„ - e<j) 



{ * + 1 , } 



dx. 



(36) 



The TTSA is defined as 
da(<p - 



A 



UT{ 



(f>s) - da(4> - 4>s + tt) 



d(7(<^ — 0s) + da{(j) — <ps + tt) 



.4 



sin(0-0 s ) 
(7T 



>)• 



(37) 



The A : 



f.P/' 



amplitude of the TTSA can be expressed 
in terms of the spin flip and spin non-flip amplitudes as 

EH: 



,sin(4>-4>s) 
1 (7T 



(38) 



Im(AB*)\A T \/m 



(i - e)\A\ 2 - (e + 



t 

4 m- 



\B\ 2 -2eRe{AB*)' 



Note that using the Trento convention [3T] the sign of this 
equation is opposite to that in Ref. |21j and the normal- 
ization is larger by a factor of tt/2. 

The cross section is calculated using both factorized 
and Regge ansatze for GPDs 5 . The value b = 1 is taken 
for the profile parameter both for sea and valence quarks. 
It is found that using b sea — oo instead of b sea = 1 leads 
to a rise of the cross section by a factor of about 1.15. The 
value b va i = b sea = 1 is chosen to provide a direct com- 
parison to previous calculations [59TI2T] . The value of the 
profile parameter for gluons is chosen as b = 2 and it has 
been checked that choosing b = 1 or oo does not change 
the cross section by more than 20%. The VF-dependence 
of the cross section for Q 2 =4 GeV 2 is shown in Fig. [7J For 
both ansatze the calculations overshoot considerably the 
experimental data from Hermes 0U|. However, a signifi- 
cant reduction of the calculated cross section might be ex- 
pected if transverse motion effects are taken into account 
|21U39| . On the other hand, also the double distribution 
based calculations of the DVCS cross section have been 
found to overshoot the data from HI |41U42j . 

An unexpected result of the calculation shown in Fig.0 
is a quite small (15-20%) pure quark contribution to the 
cross section, while in Refs. |2 11139) the quark contribu- 
tion was found to be dominant. Comparing the calcu- 
lated quark contribution to experimental data (Fig. UJ it 
could also be concluded that the gluon contribution in the 
present calculation is substantially overestimated, while 
the quark contribution itself is reasonable and can explain 
alone (in the factorized ansatz) the value of the measured 
cross section. However, there exists experimental evidence 
that the gluon spin non-flip part is indeed large EH- 

On the amplitude level, the cross sections of p° and <fi 
mesons are given as: 

<V = C p o\T q + T g \ 2 

= C p0 (\T q \ 2 + 2\T q \\T g \ cosfo> TO ) + \T g \ 2 ), (39) 
2 



The principal values of the integrals in Eqs. I35H and 13611 



are calculated in the following way: f *-^p-dx = f(b) ln(6) — 

a<0 

b 

f(a) ln(a) — J f'(x) In(a;)da: + J f'(x) ln(x)dx. In this way the 



in the subsequent calculations the exact formulae were used non-integrable singularity is exchanged by an integrable one. 
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4.2 Expected Value of TTSA and Projected Statistical 
Uncertainty 



Fig. 7. The calculated IV-dependence of the hard exclusive 
p° electroproduction cross section at Q 2 = 4 GeV 2 for the 
factorized (top) and Regge (bottom) GPD models compared 
to Hermes data |40|. T q is the quark amplitude and T g the 
gluon amplitude. The quark real part is very small and is not 
shown. 



Here T q is the quark amplitude and T g the gluon amplitude 
(the s-quark contribution to the <fr production amplitude is 
neglected) and tp qg is the effective phase between the quark 
and gluon amplitudes. In the existing GPD-based calcu- 
lations 21,32], both quark and gluon contributions are 
dominated by the imaginary parts which have the same 
sign, i.e. ip qg ~ 0. In the present calculation (p qg ~ 30° 
is obtained. Considering the wave functions of p° and <fi 
mesons to be similar (as it is supported by the measured 
values of their decay widths), C p o ~ follows, and the 
ratio of <f> to p° cross sections reads: 



0> 



\r g \ 2 



2\T q \\T g \cos(<p qg ) + \T g \ 



(40) 



At Hermes, the ratio of a^/a p o was measured |44j . 
The experimental value was 0.08±0.01, slightly increasing 
with Q 2 . Inserting it to the l.h.s of Eq. H4Q(I and taking 

<p qg = 0°(30°) yields = 0.7 (0.78). This value 

' 3 ' Hermes 

is in good agreement with the results of the present calcu- 

T I 

lation, where -pj^j ranges between 0.8 and 0.5 (0.5 and 0.3) 
for the factorized (Regge) ansatz when W increases from 
4 to 6 GeV. This is in contrast to the above mentioned 
result of a dominant quark contribution, j^j ~ 3 P51I21 - ] . 

Hence, it is concluded that H g can not be neglected, i.e. to 
arrive at the measured cross section both quark and gluon 
amplitudes have to be scaled down in a similar proportion. 



The A 



UT 



b-4>s) 



amplitude of the TTSA in Eq. (33) is 
calculated at Hermes kinematics. The statistical error 
is extrapolated from a preliminary analysis of the Her- 
mes longitudinal target-spin asymmetry measured on the 
deuteron j^SJ that is based on 8 million DIS events. In the 
latter analysis the data is not split into parts correspond- 
ing to longitudinal and transverse virtual photons, while 
the present calculation is related to longitudinal photons 
only. At Hermes kinematics {(Q 2 ) ~ 2 GeV 2 ), longitudi- 
nal photons constitute about 50% of all virtual photons. 
Also, the transverse target polarization is 0.75 while the 
longitudinal one is 0.85. The projected statistical error for 
8 million DIS events taken on a transversely polarized tar- 
get is then larger by a factor v^jppf =1.6 compared to 
that of Ref. Note that this error estimate may be 
considered 'optimistic', since it assumes that the contri- 
bution to the asymmetry from longitudinal and transverse 
photons can be completely disentangled. 



The calculated xb and ^-dependences of A 



s\n(<j>-<j>s) 
UT 



are 

shown in Fig. [H] for different values of J u . As in the case 
of DVCS, Jd is fixed inspired by the fact that the d-quark 
contribution is still suppressed, although the suppression 
in p° production is half as strong as in DVCS. Again, the 
choice of Jd = is based on the results of recent lattice 
calculations (see e.g. Ref. 36;). Note that in contrast to 
DVCS, E = results in a vanishing asymmetry. As it can 
be seen from comparing Fig. [SI to Figs. 0] and El the ex- 
pected magnitude of A s ^p~^ s ' in p° production is much 
smaller than that in DVCS. This is due to a large glu- 
onic contribution to the amplitude, which is considered 
as "passive" (E g = 0,H g ^ 0), i.e. the gluons dilute the 
asymmetry in this case. It was found that the difference 
in A B ^^ > ~^ S ' between the factorized and Regge ansatze 
is negligible. Also the variation of b sea only leads to a 
small difference as can be seen when comparing the left 
and right panels of Fig. |SJ where xb- and i-dependences 
of the A s ^^ < ^ s * > amplitude of the asymmetry are shown 
for b sea = 1 and b sea = oo, respectively. The amplitude 
of the integrated TTSA is shown in Fig. [51 for the same 
two cases. It is essentially independent of b sea and ranges 
between values of 0.10 and 0.01 when J u ranges between 
zero and 0.4. 

The projected statistical error for the integrated TTSA 
amplitudes is 0.034. Extrapolating the knowledge on the 
systematic uncertainty from Ref. |45j . its size can be ex- 
pected to be about 0.02, such that a total experimental 
uncertainty below 0.04 appears as a realistic estimate. Al- 
together, the difference in A s ^^ due to a change of 
J u between zero and 0.4 corresponds to an about 2a ef- 
fect, where a denotes the total experimental uncertainty. 
Thus it can be expected that the upcoming p° electropro- 
duction measurements performed at Hermes will provide 
an additional constraint on the size of J u . 

A tempting possibility provided by p° production is re- 
lated to an estimate of the gluonic content of E. Strongly 
simplifying, Eq. H38|l represents the ratio E/H cx (E q + 
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Fig. 8. Comparison of expected ^4^!^ ^ amplitudes of the 
p° TTSAs calculated in the Regge ansatz with b va i = 1 and 
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oo (right). Average kinematic values 



(— t) = 0.14 GeV and (xb) = 0.085 for xb and t-dependences, 
respectively, and (Q 2 ) = 2 GeV 2 correspond to a preliminary 
analysis of Hermes data on a longitudinally polarized deu- 
terium target |45|. Projected statistical errors are shown. The 
systematic uncertainty is expected to be smaller than the sta- 
tistical one. 
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Comparison of expected ^4^-^ 



amplitudes of the 



p" TTSAs calculated at average Hermes kinematics ({— t) = 
0.14 GeV 2 , (x B ) = 0.085, (Q 2 ) = 2 GeV 2 ) in the Regge ansatz 
with b va i = 1 and b sea = 1 (left) or b sea = oo (right). Pro- 
jected statistical errors are shown. The systematic uncertainty 
is expected to be smaller than the statistical one. 



Eg)/(H q + H g ). Hence, when comparing the earlier cal- 
culations |2J where gluons have been neglected (E g = 
H g = 0) to the case of "passive" gluons presented above 
[Eg = 0, Hg 7^ 0), the asymmetry gets smaller ('diluted') 
by the presence of the term containing H g in the denom- 
inator. On the other hand, if the measured asymmetry 
would be found large, this could imply that the gluons 
are "active" (E g ^ 0), so that their contribution to the 
spin-flip amplitude can not be neglected. 



5 Summary and Outlook 

Transverse target-spin asymmetries (TTSAs) in DVCS and 
p° elastic electroproduction are the only candidates known 
by now to access the GPD fiona proton target, in which 
E comes as a leading term. A code [HEl based on the model 
developed in Ref. |391T2T] is used to calculate the expected 
TTSAs to be measured in DVCS on the Hermes trans- 



versely polarized hydrogen target. To check the accessi- 
bility of E at Hermes, different parameterization ansatze 
and parameters of H and E are chosen. As the model for E 
depends on the total angular momentum of the u-quarks 
in the proton, the possibility arises to check the sensitivity 
of the data to different values chosen as J u — 0.4, 0.2, 0.0, 
while on the basis of u-quark dominance and recent lat- 
tice calculations (see e.g. Ref. |36|1 a fixed value Jd = 
is used. The calculations are performed at the Hermes 
average kinematic values [22] - The results show that the 
DVCS TTSA amplitude ^ s J»W-^) cos * is sensitive to the 
GPD E and thus to the total it-quark angular momentum 
J u , while ^4™^(* _ * s ) sln * j s n0 ^ Jt was found that aside 

from J u the amplitude A^^^^ s ' cos ^ is largely indepen- 
dent on different parameterization ansatze and model pa- 
rameters. Projected statistical errors for the asymmetries 
are evaluated by converting the ones from Ref. [22] to a 
data set corresponding to 8 million DIS events taken on a 
transversely polarized hydrogen target. 

The same parameterizations are used to calculate the 
TTSA in p° electroproduction by longitudinal virtual pho- 
tons. The main difference to the DVCS case is a large 
gluonic contribution to the amplitude. At present, only 
the spin-nonflip part of the gluonic amplitude can be rea- 
sonably described, while the spin-flip gluonic GPD E g is 
totally unknown. Therefore, throughout the calculation 
E g is set to zero ("passive" gluons). Under this assump- 
tion, the situation in p° electroproduction appears less fa- 
vorable concerning the sensitivity of the expected TTSA 
amplitude to the total angular momentum J u . However, 
should the value of the amplitude be measured larger than 
that predicted by these calculations, this would imply that 
Eg can not be neglected, and thus indicate that gluons in- 
side the proton carry significant orbital angular momen- 
tum. 

Altogether, transverse target-spin asymmetries in both 
DVCS and p° electroproduction are studied to evaluate 
projected uncertainties for extracting the value of J u from 
future data. Considering all anticipated Hermes data to 
be taken for DVCS (/9°-production), the projected total 
experimental la-uncertainty is estimated to correspond 
to a range of about 0.1 (0.2) in J u . 
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Appendix: TTSA Calculation in DVCS 

A code is used to estimate the TTSA related to DVCS. 
The coordinate system and angles defined in the code are 



10 



F. Ellinghaus et al: Can the angular momentum of u-quarks in the nucleon be accessed at Hermes? 



the same as depicted in Fig. OH The polarization of the tar- 
get in the code is defined according to the virtual photon 
direction. For a transversely polarized target, the target 
polarization direction can be chosen either in the lepton 
plane (x direction) or perpendicular to it (y direction). 
The former corresponds to <fis = or tt, the latter to 
4>s = tt/2 or 3tt/2. Therefore the following intermediate 
asymmetries can be calculated: 



A x (4>) 



Ay(4>) 



da,p s= o(<j)) - da < p s= ^{4)) 
da,p s=0 ((j}) + da,p s=7T (<j)) ' 
do> s =f {4>) - da, 3^ (0) 



duA 



-WO 



(<py 



Defining the following functions 

A^cj)) = A x ■ smcj)- 
A 2 (4>) = A x ■ cos^- 



A y ■ COS ( 

- A, ■ sin < 



(41) 



(42) 



the contribution of the transverse target polarization com- 
ponent of the interference term Ttp to the total cross sec- 
tion in Eq. (|25D can be expressed as: 



dai 



da,. 



Ax (</>)• sin (<f> - <j) S )+ A 2 ((f>) ■ cos (0 - </> s ) 



(43) 

Therefore the asymmetries defined in Eq. (|28() can be com- 
puted as: 



sin (4>—(fis) cos fp 
UT 

cos (4>—'fis) s i n 4 1 
UT 



A COS (f) 
^1 5 

a sin (p 
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